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Abstract 



We consider the Ginzburg-Landau Hamiltonian with a cubic-symmetric quar- 
tic interaction and compute the renormalization-group functions to six-loop 
order in d = 3. We analyze the stability of the fixed points using a Borel 
transformation and a conformal mapping that takes into account the singu- 
larities of the Borel transform. We find that the cubic fixed point is stable for 
N > N c , N c = 2.89(4). Therefore, the critical properties of cubic ferromag- 
nets are not described by the Heisenberg isotropic Hamiltonian, but instead by 
the cubic model at the cubic fixed point. For N = 3, the critical exponents at 
the cubic and symmetric fixed points differ very little (less than the precision 
of our results, which is < 1% in the case of 7 and v). Moreover, the irrelevant 
interaction bringing from the symmetric to the cubic fixed point gives rise to 
slowly-decaying scaling corrections with exponent 0J2 = 0.010(4). For N = 2, 
the isotropic fixed point is stable and the cubic interaction induces scaling 
corrections with exponent L02 = 0.103(8). These conclusions are confirmed by 
a similar analysis of the five-loop e-expansion. A constrained analysis which 
takes into account that N c = 2 in two dimensions gives N c = 2.87(5). 

PACS Numbers: 75.10.Hk, 05.70.Jk, 64.60.Fr, ll.10.Kk 
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I. INTRODUCTION 



According to the universality hypothesis, critical phenomena can be described in terms 
of quantities that do not depend on the microscopic details of the system, but only on 
global properties such as the dimensionality and the symmetry of the order parameter, and 
the range of the interactions. There exist several physical systems that are characterized 
by short-range interactions and an iV-component order parameter. Because of universality, 
their critical properties can be studied by using the Ginzburg-Landau 4 Hamiltonian and 
by employing standard field-theoretic renormalization-group techniques. When the order 
parameter has only one component, one obtains the Ising universality class that describes 
for instance the liquid-vapor transition in simple fluids and the transitions of multicom- 
ponent fluid systems; in this case the density plays the role of the order parameter. The 
two-component model (XY model) describes the helium superfluid transition, the Meissner 
transition in type-II superconductors and some transitions in liquid crystals, while the limit 
iV — > gives the infinite- length properties of dilute polymers in a good solvent. 

The critical properties of many magnetic materials are also computed using the N- 
component Ginzburg-Landau Hamiltonian. Uniaxial (anti-)ferromagnets should be de- 
scribed by the Ising universality class (N = 1), while magnets with easy-plane anisotropy 
should belong to the XY universality class. Ferromagnets with cubic symmetry are often 
described in terms of the N = 3 Hamiltonian. However, this is correct if the non-rotationally 
invariant interactions that have only the reduced symmetry of the lattice are irrelevant in the 
renormalization-group sense. Standard considerations based on the canonical dimensions of 
the operators indicate that there are two terms that one may add to the Hamiltonian and 
that are cubic invariant: a cubic hopping term Z^=i,3(c^</v) 2 an d a cubic interaction term 
J2ix=i,3 <f>t,- The first interaction was studied in Refs. |l]-|5|. A two-loop 0(e 2 ) calculation 
indicates that it is irrelevant at the symmetric point, although it induces slowly-decaying 
crossover effects. We will not consider it here — although it would be worthwhile to perform 
a more systematic study — since the second term already introduces significant changes in 
the critical behavior of the system. We will therefore consider a three-dimensional 4 theory 
with two quartic couplings 



The added cubic term breaks explicitly the O(N) invariance of the model, leaving a residual 
discrete cubic symmetry given by the reflections and permutations of the field components. 

The model described by the Hamiltonian has been extensively studied. It has four 
fixed points [|6|||: the trivial Gaussian one, the Ising one in which the N components of the 
field decouple, the 0(iV)-symmetric and the cubic fixed points. 

The Gaussian fixed point is always unstable, and so is the Ising fixed point [[/]]. Indeed, in 
the latter case, it is natural to interpret Eq. ( |1 . 1| ) as the Hamiltonian of N Ising-like systems 
coupled by the 0(iV)-symmetric term. But this interaction is the sum of the products of 
the energy operators of the different Ising systems. Therefore, at the Ising fixed point, 
the crossover exponent associated to the 0(iV)-symmetric quartic term should be given by 
the specific-heat critical exponent ai of the Ising model, independently of N. Since ai is 
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FIG. 1. Renormalization-group flow in the coupling plane (u, v) for N < N c and N > N c . 

positive, indeed oti = 0.1099(7) (see e.g. Ref. |J and references therein), the Ising fixed 
point is unstable. 

While the Gaussian and the Ising fixed points are unstable for any number of components 
N, the stability properties of the 0(iV)-symmetric and of the cubic fixed points depend on 
N. For sufficiently small values of N, N < N c , the 0(iV)-symmetric fixed point is stable and 
the cubic one is unstable. For iV > N c , the opposite is true: the renormalization-group flow 
is driven towards the cubic fixed point, which now describes the generic critical behavior 
of the system. The 0(iV)-symmetric point corresponds to a tricritical transition. Figure [I] 
sketches the flow diagram in the two cases N < N c and N > N c . At N = N c , the two 
fixed points should coincide, and logarithmic corrections to the 0(iV)-symmetric critical 
exponents are expected. Outside the attraction domain of the fixed points, the flow goes 
away towards more negative values of u and/or v and finally reaches the region where the 
quartic interaction no longer satisfies the stability condition. These trajectories should be 
related to first-order phase transitions P,[IT|. 

If TV > N c , the cubic anisotropy is relevant and therefore the critical behavior of the 
system is not described by the Heisenberg isotropic Hamiltonian. If the cubic interaction 
favors the alignment of the spins along the diagonals of the cube, i.e. for a positive coupling 
vq, the critical behavior is controlled by the cubic fixed point and the cubic symmetry is 
retained even at the critical point. On the other hand, if the system tends to magnetize along 
the cubic axes — this corresponds to a negative coupling t> — then the system undergoes a 
first-order phase transition |T5|j5|, |l^ , p^1 . Moreover, since the symmetry is discrete, there are 
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no Goldstone excitations in the low-temperature phase. The longitudinal and the transverse 
susceptibilities are finite for T < T c and H — > 0, and diverge as |t| -7 for t oc T — T c — > 

In the limit N — > oo, keeping iVtt and t> fixed, one can derive exact expressions for the 
exponents at the cubic fixed point. Indeed, in this limit the model can be reinterpreted 
as a constrained Ising model [pOR , leading to a Fisher renormalization of the Ising critical 
exponents One has 
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where rji, Vi, and are the critical exponents of the Ising model (see e.g. 
references therein for recent estimates of the Ising critical exponents). 

If N < N c , the cubic term in the Hamiltonian is irrelevant, and therefore, it generates 
only scaling corrections |t| Ac with A c > 0. However, their presence leads to important 
physical consequences. For instance, the transverse susceptibility at the coexistence curve 
(i.e. for T < T c and H — > 0), which is divergent in the 0(iV)-symmetric case, is now finite 



and diverges only at T c as |t| _7_Ac [^,^,^,|,^]. In other words, below T c , the cubic term 
is a "dangerous" irrelevant operator. Note that for iV sufficiently close to N c , irrespective 
of which fixed point is the stable one, the irrelevant interaction bringing from the unstable 
to the stable fixed point gives rise to very slowly decaying corrections to the leading scaling 
behavior. 

In three dimensions, a simple argument based on the symmetry of the two-component 
cubic model |26| shows that the cubic fixed point is unstable for N = 2, so that N c > 2. 
Indeed, for N = 2, a ir/4 internal rotation, i.e. 



+ 



-02), 



;i.3) 



maps the cubic Hamiltonian ( |1.1| ) into a new one of the same form but with new couplings 

KX) § iven b y 



u + § v , 



-v . 



(1.4) 



This symmetry maps the Ising fixed point onto the cubic one. Therefore, the two fixed points 
describe the same theory and have the same stability. Since the Ising point is unstable, the 
cubic point is unstable too, so that the stable point is the isotropic one. In two dimensions, 
this is no longer true. Indeed, one expects the cubic interaction to be truly marginal for 
N = 2 f27]j28[l and relevant for iV > 2 [29|, so that N c = 2 in two dimensions. 

During the years, the model (|1 . 1|) has been the object of several studies 
P J2^|j3T| -p3|,p|,^,|34]- p^ , p^ ,p^-p6| . In the 70's several computations were done using the e- 
expansion |6| j23| , |32|j33|] ; they predicted 3 < N c < 4, indicating that cubic ferromagnets are 
described by the 0(iV)-invariant Heisenberg model. However, recent studies have ques- 
tioned these conclusions. Field-theoretic studies, based on the analysis of the three-loop 



3T| , |3"8| and four-loop series |I9],[46| in fixed dimension, and of the five-loop expansion in pow- 
ers of e = 4 — d p0|- |43|j46[| suggest that N c < 3, although they do not seem to be conclusive 
in excluding the value iV c = 3. On the other hand, the results of Ref. [|44||, obtained from 
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TABLE I. Summary of the results in the literature. The values of the exponents refer to 
N = 3. The subscripts "s" and "c" indicate that the exponent is related to the symmetric and to 
the cubic fixed point respectively. 



method 



results 



Ref. |2 



Ref. p4 



Ref. |6 



Ref. [28 



Ref. |9 



Ref. gO 



Ref. (42 



Ref. g3 



Ref. 144 



Ref. |6 



1974 
1977 
1981 
1982 
1989 
1995 
1997 
1997 
1998 
1999 



This work 
This work 



e-exp.: 0(e 3 ) 
approx. RG 

H.T. exp.: 0(/3 10 ) 
scaling-field 

d = 3 exp.: 0(g 4 ) 
e-exp.: 0(e 5 ) 
e-exp.: 0(e 5 ) 
e-exp.: 0(e 5 ) 
Monte Carlo 

d = 3 exp.: 0(g A ) 
e-exp.: 0(e 5 ) 

d = 3 exp.: C% 6 ) 



iV c ~ 3.128 
i/,W2, a = -0.11, iV c ~2.3 
i/ s W2,. = -0.63(10), iV c <3 
A c ~ 3.38 
lo 2 ,c ^ 0.008, N c ~ 2.91 
iV c ~ 2.958 
u%s = -0.00214, uj 2fi = 0.00213, N c < 3 
N c ~ 2.86 
w 2 ,s = 0.0007(29), N c « 3 
W2ifl = -0.0081, w 2iC = 0.0077, A c = 2.89(2) 
W2)S = -0.003(4), u; 2 , c = 0.006(4), A c = 2.87(5) 
W2)S = -0.013(6), w 2 , c = 0.010(4), A c = 2.89(4) 



Monte Carlo simulations using finite-size scaling techniques, are perfectly consistent with 
the value N c ~ 3. The authors of Ref. [P] even suggest that iV c = 3 exactly. 

A further study of this issue is therefore of particular relevance for the ferromagnetic ma- 
terials characterized by an order parameter with = 3. For this purpose we extended the 
perturbative expansions of the /3-functions and of the exponents to six loops in the frame- 



work of the fixed-dimension field-theoretic approach |47| • These perturbative expansions are 
only asymptotic. Nonetheless, accurate results can be obtained by employing resummation 
techniques that use their Borel summability |48[] and the knowledge of the large-order be- 
havior [4^,50]. For this reason, we have also computed the singularity of the Borel transform 
that is closest to the origin, extending the calculations of Refs. [ff9j , [50[| . 

The analysis of the perturbative series has been done following closely Ref. [^J. We have 
estimated errors using an algorithmic procedure, trying to be as conservative as possible. 
This can be immediately realized by comparing our uncertainties with those previously 
quoted: even though our series are longer, the errors we report are sometimes larger than 
those of previous studies. Our results confirm previous field-theoretic studies: the A^ = 3 
isotropic fixed point is indeed unstable and we estimate A^ = 2.89(4) from the six- loop fixed- 
dimension expansion and N c = 2.87(5) from the reanalysis of the five-loop e-expansion. For 
comparison, in Table [I] we report our estimates together with previous determinations of 
Ac and of the eigenvalues for A^ = 3. It should be noted that the estimates of the critical 
exponents do not essentially depend on which fixed point is the stable one. Moreover, the 
tiny difference (smaller than the precision of our results, which is ^ 1% in the case of 7 and 
v) between the values at the two fixed points would be very difficult to observe, because of 
crossover effects decaying as t A with A = c^.c^c = 0.007(3). Large corrections to scaling 
appear also for N = 2. Indeed, at the XY fixed point (the stable one), we find u 2 = 0.103(8). 
Thus, even though the cubic interaction is irrelevant, it induces strong scaling corrections 
behaving as t A , A = uj 2 v ~ 0.06. Therefore, crossover effects are expected in this case, 
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depending on the strength of the cubic term. Finally, we have checked the theoretical 
predictions for the model in the large- N limit finding good agreement. 

We want to mention that, in the limit iV — > 0, the cubic model ( |1 . 1| ) describes the Ising 



model with site-diluted disorder |p2|-j54|j . However, in this case, the perturbative expansion is 
not Borel summable ]55|-p7| . Therefore, it is not completely clear how to obtain meaningful 
results from the perturbative series. An investigation of these problems will be presented 
elsewhere. 

The paper is organized as follows. In Sec. |I|we present our calculation of the perturbative 
expansions to six loops in d — 3. We give the basic definitions, the six-loop series, and the 
singularity of the Borel transform. In Sec. [TTI| we present the analysis of these expansions: 
we determine the stability of the fixed points and compute the exponents for several values 
of N. In Sec. [TV] we present a reanalysis of the e-expansion five-loop series. The new 
analysis differs from the previous ones in the fact that it uses the large-order behavior of the 
series at the cubic fixed point. Finally, in the Appendix we report a three- loop e-expansion 
computation of the zero-momentum four-point couplings at the cubic fixed point in three 
dimensions. 



II. THE FIXED-DIMENSION PERTURBATIVE EXPANSION IN THREE 

DIMENSIONS 

A. Renormalization of the 4 theory with cubic anisotropy 

The fixed-dimension 4 field-theoretic approach j|7|] provides an accurate description of 



the critical properties of 0(iV)-symmetric models in the high-temperature phase (see e.g. 
Ref. H58| ). The method can also be applied to the two-parameter cubic model [^7f. The 



idea is to perform an expansion in powers of appropriately defined zero-momentum quartic 
couplings. In order to obtain estimates of the universal critical quantities, the perturbative 
series are resummed exploiting their Borel summability, and then evaluated at the fixed- 
point values of the couplings. 

The theory is renormalized by introducing a set of zero-momentum conditions for the 
(one-particle irreducible) two-point and four-point correlation functions: 

rS(p) = s^ 1 U 2 + P 2 + o(p A )} , (2.i) 



r2L(o) = V™ 



u 



(Sabred + 5 ac 5 bd + 5 ad 5 bc ) + v 8 ab 8 ac 8 ad 



(2.2) 



They relate the second-moment mass m, and the zero-momentum quartic couplings u and 
v to the corresponding Hamiltonian parameters r, uq and Vq: 



u Q = muZuZ^ , v = mvZ v Z (j) . 

In addition, one introduces the function Z t that is defined by the relation 



(2.3) 



r a 1 b 2 \o) = s ab z ) 



-1 
t ' 



(2.4) 
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where I^ 1,2 ) is the (one-particle irreducible) two-point function with an insertion of \4> 2 - 
From the pertubative expansion of the correlation functions T^ 2 \ and I^ 1 ' 2 ) and the 

above relations, one derives the functions Z<f,(u,v), Z u (u,v), Z v (u,v), Z t (u,v) as a double 

expansion in u and v. 

The fixed points of the theory are given by the common zeros of the /^-functions 

(J u {u,v) = m — 
am 



(J v {u,v) = m — 
am 



U ,VQ 



U ,VQ 



(2.5) 



The stability properties of the fixed points are controlled by the eigenvalues Ui of the matrix 



( d/3 u (u,v) d/3 u (u,v) \ 

du dv \ 

df3 v (u,v) d(3 v {u,v) I ' 

du dv ' 



(2.6) 



computed at the given fixed point: a fixed point is stable if both eigenvalues are positive. 
The eigenvalues uii are related to the leading scaling corrections, which vanish as 
where Aj = vuii. 

One also introduces the functions 



T] t (u,v) 



d\nZ, 



dlnm 
dh\Z t 



d\\im 



u ,v 



u {) .ro 



dlnZ^ dlnZ, 

" Pu T Pv~ 

ou 

Pu 71 fiv 



dv 
d\nZ t 



du 



dv 



Finally, the critical exponents are obtained from 

V = V<t>( u *, v *), 

u=[2- V<i> (u*,v*)+m(u*,v*)]- 1 
7 = 1/(2-77). 



■ fI 

(2.7) 
(2.8) 



(2.9) 
(2.10) 
(2.11) 



B. The six-loop perturbative series 



We have computed the perturbative expansion of the correlation functions (j2.1p, ( [2.2| ) 
and (p.4j) to six loops. The diagrams contributing to the two-point and four-point functions 
to six-loop order are reported in Ref. ||59|| : they are approximately one thousand, and it is 
therefore necessary to handle them with a symbolic manipulation program. For this purpose, 
we wrote a package in Mathematica It generates the diagrams using the algorithm 
described in Ref. [6~I], and computes the symmetry and group factors of each of them. 



We did not calculate the integrals associated to each diagram, but we used the numerical 
results compiled in Ref. p9f . Summing all contributions we determined the renormalization 



constants and all renormalization-group functions. 

We report our results in terms of the rescaled couplings [B2 



7 



167T _ 167T _ 

u = — R N u, v = — v, (2.12) 

where iiV = 9/(8 + N), so that the /3-functions associated to u and v have the form 
(3-(u, 0) = -u + u 2 + 0(u 3 ) and /3 e (0, v) = -v + v 2 + 0(w 3 ). 
The resulting series are: 

, 2 4(190 + 41AH , 

p a = -u + u 2 + -uv- --— ttt^m 3 (2.13) 



3 27(8 + N) 



400 2 92 __ 2 (u) 

; tU V UV +U > 0)a 

81(8 + N) 729 ^ " 



i+j>3 

2 12 308 o 832 2 , 

/3 = -v + V 2 H — ?T ; T-rrUV (2.14) 

M 8 + N 729 81(8 + N) v ' 



4(370 + 23AQ _ 2 _ _ _ 

; r-r— tt V + f > /uV, 

27(8 + AT) 2 « 



8(2 + AH „ 16 8 , ^ ™ , , , 

"* = 2f{sTW a + 8l(8T7V) fiB + 729 s + J>« W ' (2 ' 15) 

(2 + AH 1 2(2 + AT) „ 4 2 „ „ m ■ , , 

* = -JsTWf ~ f> + (W" + WTn) m + 27 11 + J>8 W - < 2 - 16 > 

For 3 < z + j < 6, the coefficients 6^, b[j\ ej^ and ej^ are reported in Tables [TJ, [TIT], |IV] 
and respectively. 

We have performed several checks of our calculations: 

(i) Pu(u,0), T](j, (u, 0) and ^(w, 0) reproduce the corresponding functions of the O(N)- 
symmetric model []6~2|j6~3l ; 

(ii) /3y(0, v), r),j,(0,v) and rjt(0,v) reproduce the corresponding functions of the Ising-like 
(N = 1) 4 theory; 

(iii) The following relations hold for N = 1: 

(3u{u,x -u) + /3a(u,x - u) = f3y(0,x), (2.17) 
r^(w,x - u) = 77^(0, x), 
r) t (u,x -u) = T}t(0,x). 



(iv) For N = 2, using the symmetry Q1.3Q and ( |1.4|) , and taking into account the rescalings 
( |2.12j ), one can easily obtain the identities 



P a (u + §v, -v) + + §«, -v) = ^(w, u), (2.18) 



5 
3 

/3 s (u + |u, -u) = u), 
r/t(iZ + |u, -«) = Vt(u,v). 
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TABLE II. The coefficients b\y , cf. Eq. ( gig ). 



hi 


R7/f> (u) 

IV 27 




3 


(1 97*38^17 4- 07^61090 V -4- 0018^01016 A/" 2 
u. z i oojiji i t^ u.ui jou'iu^y i v -p u . uuioju'iuiu i v 




2 1 


6774939^ -U 0973^400 /V 
u.ui i^iZrO^u t^ u.uz i ojoiuy i v 




1.2 


41 5456^ 4- 002 c i l SQ21 48 /V 




3 


0Q0448QM 




4 o 


970,9^794 0Q1 8'3 4: i74Q AT 00 "VI ^0^6 -1 6 N' A 4- fl 0000937998Q'} /V d 
— u.z i yzo ( — u.uyiooo/ i iy iv — u.uuu'iuy uu^u i v -\- u.uuuuzo i zzoyo i v 




3 1 


—0 94S8S662 — 08S252807 TV 4- 00061860174 A 7 " 2 




2 2 


—0 Q64Q7888 — 01946014^ /V 




13 


—0 42331 874 — 001 770Q42Q N 




4 


—0 07^4466Q9 

— U.UI ij'± £ ±UUyZ 






^^17/1477-1-0 1 3949 ^09 A/" -L 01 1 '399098 AT'^ -L 0000^ W^'^71 V'" 5 ' si fi7f : ;SQQQ 1 0" 

u.ooi / ■rt^t Mt" U.10/4ZOUZ iv -r u.ui iozzuzo iv T" u.uuuuo^tooo / ±y i v T"O.Di uoyoo 1U 


IV 


i i 


1 ^900008 j_ fl 1 Q4^0^^fi AT _L 001 1 07881 4 Af2 i n 0000^1 77Q789 Af3 

i .ozuyuuo T" u. ly^touooD iv t* u.uujliu i oui^t i v -t- u.uuuuoi 1 1 y i oz i v 






9 9073'347 _l_ Ofi c l'3'ifi'}98 A/" -L 000^^8409^ Af2 
Z.ZU 1 Oo4 l T" U.UO DOOUOZO iV j U.UUUO 0044 y ZO iv 




2,3 


1.5315693 + 0.010676901 N 




1,4 


0.56035196 + 0.0013469481 N 




0,5 


0.087493302 




6,0 


-0.51049889 - 0.21485252 N - 0.023839375 - 0.00050021682 N a + 2.0167763 10 


-o N* + 4.4076733 10" s N° 


5,1 


-2.6984083 - 0.45068252 N - 0.010821468 N 2 + 0.00005796668 TV 3 + 2.0515456 10~ 




4,2 


-5.1135549 - 0.26769177 N - 0.0006311751 N 2 + 0.000019413374 N 3 




3.3 


-4.9317312 - 0.067574712 N + 0.000028278087 N 2 




2,4 


-2.754683 - 0.0095836704 N 




1,5 


-0.86229463 - 0.001856332 N 




0,6 


-0.1179508 





These relations are exactly satisfied by our six- loop series. Note that, since the Ising fixed 
point is (0, <?!) with g\ = 1.402(2) ||, the above symmetry gives us the location of the cubic 
fixed point: (§#J,-#|). 

(v) In the large-iV limit the critical exponents of the cubic fixed point are related to 
those of the Ising model: r\ = r/j and u = uj / (1 — «/). One can easily see that, for iV — > oo, 
rj^UjV) = rji(v), where 777(f) is the perturbative series that determines the exponent r\ of 
the Ising model. Therefore, the first relation is trivially true. On the other hand, the second 
relation v = — oti) is not identically satisfied by the series, and is verified only at the 

critical point |M . 



We finally note that our series are in agreement with the five-loop results appeared 



recently in Ref. [55 



C. Singularity of the Borel transform 

Since field-theoretic perturbative expansions are asymptotic, the resummation of the 
series is essential to obtain accurate estimates of the physical quantities. In three-dimensional 
</> 4 theories one exploits their Borel summability |58| and the knowledge of the large-order 
behavior of the expansion (see e.g. Ref. fl58|1 ). 

In the case of the 0(iV)-symmetric </> 4 theory, the expansion is performed in powers of the 
zero-momentum four-point coupling g. The large-order behavior of the series S(g) = J2 Sk9 k 
of any quantity is related to the singularity g>& of the Borel transform that is closest to the 
origin. Indeed, for large k, 

s k ~ k\ (~a) k k b \l + Oik" 1 )] with a = -l/g b . (2.19) 
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TABLE III. The coefficients b\f , cf. Eq. 



h j 




N 2.7 


3,0 


0.64380517 


+ 0.05741276 TV - 0.0017161966 TV 2 


2,1 


1.6853305 4 


0.0030714114 TV 


1,2 


1.3138294 




0,3 


0.3510696 





4~0 -0.76706177 - 0.089054667 TV + 0.000040711369 W - 0.000087586118 N 3 

3.1 -2.7385841 - 0.049218875 TV - 0.00002623469 TV 2 

2.2 -3.3477204 + 0.0075418394 TV 

1.3 -1.8071874 
0,4 -0.37652683 



5fi 1.0965348 + 0.15791293 TV + 0.0023584631 N' z - 0.000061471346 TV 3 - 5.3871247 10~ u TV 4 

4.1 4.9865485 + 0.17572792 TV - 0.0020718369 TV 2 - 0.000019382912 TV 3 

3.2 8.3645284 + 0.0039620562 TV + 0.00021363122 N 2 

2.3 6.8946012 - 0.0230874 TV 

1.4 2.8857918 

(V5 0.49554751 

67) -1.7745533 - 0.30404316 TV - 0.0094338079 N 2 + 0.000066993864 N 6 - 6.5724895 10"° TV 4 - 3.753114 10" ' N" 

5.1 -9.8298296 - 0.53384955 N + 0.0022033252 TV 2 - 0.00013066822 TV 3 - 2.5959429 10" 6 TV 4 

4.2 -21.073538 - 0.16628697 N — 0.000014827682 N 2 + 4.4988524 10" 6 N z 

3.3 -23.569724 + 0.095716867 N - 0.00083903999 N 2 

2.4 -14.927998 + 0.0486813 N 

1.5 -5.1298717 
0,6 -0.74968893 



TABLE IV. The coefficients e>J J , cf. Eq.(^T5|). 



i,j RljeW 

• J Z.7 

3fi 0.00054176134 + 0.00033860084 N + 0.000033860084 W 

2.1 0.002437926 + 0.00030474076 N 

1.2 0.0027426668 

0,3 0.00091422227 

4~0 0.00099254838 + 0.00070251807 N + 0.0001018116 N 2 - 6.5516886 10" ' N A 

3.1 0.0059552903 + 0.0012374633 N - 7.8620264 10" 6 N 2 

2.2 0.01046567 + 0.00031166693 2V 

1.3 0.0071848915 

0,4 0.0017962229 

5fi -0.00036659735 - 0.0002572117 TV — 0.000032026611 N 2 + 2.2430702 lO" 13 N' 6 - 1.1094045 10"' TV 4 

4.1 -0.0027494801 - 0.00055434769 TV + 0.000036974267 TV 2 - 1.6641067 10~ 6 TV 3 

3.2 -0.0064696069 - 0.000070210701 TV + 2.7823882 10 -6 TV 2 

2.3 -0.0066046286 + 0.00006759333 TV 

1.4 -0.0032685176 

0,5 -0.00065370353 

6fi 0.00069568037 + 0.00056585941 TV + 0.00012057302 TV 2 + 5.7466979 10~ B TV 3 - 3.8385183 10~ s TV 4 - 1.0441273 10" s TV b 

5.1 0.0062611234 + 0.001962173 TV + 0.00010407066 TV 2 - 3.1504745 10~ 7 TV 3 - 1.8794292 10" 7 TV 4 

4.2 0.018957129 + 0.0018483045 TV + 0.000012482494 TV 2 - 7.5537784 10" 7 TV 3 

3.3 0.027158805 + 0.00059570458 TV + 1.7043408 10~ 6 TV 2 

2.4 0.020775681 + 0.000041479576 TV 

1.5 0.0083268641 
0,6 0.0013878107 
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TABLE V. The coefficients ef] , cf. Eq. ( gig ) 



*>i 




JV 2.7 


3,0 


-0.025120499 


- 0.016979919 N - 0.0022098349 


2,1 


-0.11304225 - 


- 0.019888514 N 


1.2 


-0.13037154 - 


- 0.0025592148 N 


0,3 


-0.044310253 





3,1 
2.2 

1,3 
0,4 



0.12876028 4 
0.22779178 + 
0.15630733 + 
0.039519569 



0.029764853 N - 
0.0093256378 N 
0.0017709429 N 



0.00044686275 N 2 



- 0.00013566164 
0.000025489635 N 3 



1.699309 10" b Af 4 



5,0 
4,1 
3.2 
2,3 
1,4 
0.5 



-0.022694287 - 0.017985168 N - 0.0035835384 N 2 
-0.17020715 - 0.049785186 N - 0.0019839454 N 2 - 
-0.40917573 - 0.034538379 N - 0.00028943185 N 2 
-0.43464785 - 0.0093557007 N 
-0.22065482 - 0.0013469481 N 
-0.044400355 



0.029450619 + 0.024874579 N + 0.005728397 N 2 + 0.00031557863 N a - 

0.26505557 + 0.091343427 N + 0.0058838593 N 2 - 0.00010172194 N 3 - 

0.80694662 + 0.09798244 N + 0.00053192615 N 2 - 0.000012819748 N 3 

1.1638111 + 0.043471034 N - 0.000017867101 N 2 

0.89481393 + 0.010634231 N 

0.36032293 + 0.001856332 N 

0.060363211 



5.858689 10"° Af 4 - 
1.8672312 10~ 6 AT 4 



1.0373506 10"' N b 



6,0 
5,1 
4,2 
3,3 
2,4 
1,5 
0,6 



The value of depends only on the Hamiltonian, while the exponent b depends on which 
Green's function is considered. If the perturbative expansion is Borel summable, then is 
negative. The value of can be obtained from a steepest-descent calculation in which the 
relevant saddle point is a finite-energy solution (instanton) of the classical field equations 
with negative coupling P^J5C| . If the Borel transform is singular for g = g>&, its expansion 
in powers of g converges only for \g\ < \gjy\. An analytic extension can be obtained by a 



conformal mapping [pl|, such as 



> — g/gt— i 

y(g) = \ — • (2.20) 
-g/9b + i 



In this way the Borel transform becomes a series in powers of y(g) that converges for all 
positive values of g provided that all singularities of the Borel transform are on the real 
negative axis [pT|| . For the 0(A^)-symmetric theory accurate estimates (see e.g. Ref. pS| ) 
have been obtained resumming the available series: the /3-function ]62|] is known up to six 
loops, while the functions 77^ and rj t are known to seven loops [|S7[| . A subtle point in this 



method is the estimate of the uncertainty of the results. Indeed the non-analyticity of the 
Callan-Symanzik /5-function at the fixed-point value g* |47 , |68 -f70H may cause a slow conver- 



gence of the estimates to the correct fixed-point value. This may lead to an underestimate 
of the uncertainty which is usually derived from stability criteria. The reason is that this re- 
summation method approximates the /3-function in the interval [0, g*} with a sum of analytic 
functions. Since, for g = g*, the /3-function is not analytic, the convergence at the endpoint 
of the interval is slow. However, the comparison of these results with those obtained in other 
approaches shows that the above nonanalyticity causes only very small effects, which are 
negligible in most cases. See Refs. [170,0] for a discussion of this issue. 
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In order to apply a resummation technique similar to that used in Ref. |5T[ to our six- 
loop series, i.e. in order to use a Borel transformation and a conformal mapping to get a 
convergent sequence of approximations, we extended the large-order analysis to the cubic 
model. In particular we considered the double expansion in u and v at fixed z = v/u. Then, 
we studied the large-order behavior of the resulting expansion in powers of u. This was done 



following the standard approach described for example in Refs. p0| , |58| , i.e. by studying the 
saddle-point solutions of the cubic model. We will report the calculation elsewhere |72|| ; here 
we give only the results. 

For z = v/u fixed and N > 1, the singularity of the Borel transform closest to the origin, 
Ub, is given by 

1 -a(R N + z) for < z, (2.21) 



2NR 



±- = -a(R N + ±z) for 0>z> ^ 



where 



a = 0.14777422..., R N = . (2.22) 

N + 8 K J 

Note that the series in powers of u keeping z fixed is not Borel summable for u > and 
z < — Rn • This fact will not be a real limitation for us, since we will only consider values of 
z such that tit, < 0. It should be noted that these results do not apply to the case N = 0. 
Indeed, in this case, there are additional singularities in the Borel transform |55-57[ 



The exponent b in Eq. ( p.!9| ) is related to the number of symmetries broken by the 
classical solution [^(J. It depends on the quantity considered. In the cubic model, for v ^ 0, 
we have b = 5/2 for the function rj^, and b = 7/2 for the /3-functions and rj t . For v — 0, we 
recover the results of the 0(N) -symmetric model, that is b = 2 + N/2 for rj^, and b = 3 + N/2 
for the /3-function and r\ t |5l| . 



III. ANALYSIS OF THE FIXED-DIMENSION SIX-LOOP SERIES 

In this Section we present the analysis of the six-loop series in fixed dimension. The 
resummation of the series has been done using the method of Ref. |H] and the expression 
for the Borel singularity we have given in the previous Section. Explicitly, given a series 

R(u,v)= ^2J2 R hkU h v k , (3.1) 

k=0h=0 

we have generated a new quantity E(R) p (a,b;u,v) according to 

E(R) p (a,b;%v) = £ B k (a,b;v/u) f°° dtt b e^ (3.2) 

where 
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TABLE VI. Estimates of the fixed-point value u* and of the critical exponent UJ2 at the 
0(A)-symmetric fixed point for N = 3. The number p indicates the number of loops included in 
the analysis, the columns labelled by a and b indicate the intervals of a and b used, "final" reports 
our final estimate for the given analysis. 



a 


b 


p = 4 


p = 5 


p = 6 


final 


u* [0.5,2.5] 


[5,11] 


1.400(6) 


1.392(4) 


1.392(1) 


1.392(2) 


[0.5,2.5] 


[5,13] 


1.400(6) 


1.392(4) 


1.392(1) 


1.392(2) 


[0.5,2.5] 


[3,13] 


1.406(17) 


1.389(7) 


1.394(4) 


1.394(8) 


[-0.5,3.5] 


[5,11] 


1.404(13) 


1.392(6) 


1.393(2) 


1.393(4) 


cj 2 [0.5,2.5] 


[5,11] 


-0.014(10) 


-0.009(5) 


-0.013(2) 


-0.013(4) 


[0.5,2.5] 


[5,13] 


-0.012(9) 


-0.009(4) 


-0.012(2) 


-0.012(4) 


[0.5,2.5] 


[3,13] 


-0.017(22) 


-0.007(12) 


-0.014(7) 


-0.014(14) 


[-0.5,3.5] 


[5,11] 


-0.015(16) 


-0.008(7) 


-0.013(3) 


-0.013(6) 



TABLE VII. Estimates of u* and of the eigenvalue io<i at the 0(A)-fixed point for several 
values of N and for different orders p of the perturbative series. The last column reports the final 
estimate. 



N 


a 


b 


p = 4 


p = 5 


p = 6 


final 


— * 

u 


2 


[-0.5,3.5] 


[5,11] 


1.422(15) 


1.409(7) 


1.408(2) 


1.408(4) 


3 


[-0.5,3.5] 


[5,11] 


1.404(13) 


1.392(6) 


1.393(2) 


1.393(4) 


4 


[0.5,4.5] 


[9,15] 


1.372(12) 


1.375(2) 


1.375(1) 


1.375(2) 


8 


[0.5,4.5] 


[9,15] 


1.303(6) 


1.304(1) 


1.305(1) 


1.305(2) 



CO 1 



0J 2 

~2 [-0.5,3.5] [5,11] 0.099(20) 0.107(11) 0.103(4) 0.103(8) 

3 [-0.5,3.5] [5,11] -0.015(16) -0.008(7) -0.013(3) -0.013(6) 

4 [0.5,4.5] [9,15] -0.105(10) -0.109(3) -0.111(2) -0.111(4) 
8 [0.5,4.5] [9,15] -0.371(8) -0.379(4) -0.385(4) -0.385(8) 
co —1 



— x/u b (z) — 1 

y(x; z) = === , (3.3) 

yjl ~ X/U b (z) + 1 

and Ub(z) is denned in Eq. ( |2.21|) . The coefficients B k (a,b;v/u) are determined by the 
requirement that the expansion of E(R) p (a,b;u,v) in powers of u and v gives R(u,v) to 
order p. For each value of a, b, and p, E(R) p (a, b;u,v) provides an estimate of R(u, v). 

First of all, we have analyzed the stability properties of the 0(V)-symmetric fixed point. 
Since d{3$/du(u, 0) = 0, the eigenvalues are simply 

^1=^(^,0), u 2 =^(u*,0), (3.4) 

where u* is the fixed-point value of u. The exponent lo\ is the usual exponent that is 
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considered in the 0(iV)-symmetric theory, while u 2 is the eigenvalue that determines the 
stability of the fixed point. 

In order to compute u> 2 , for many choices of the four parameters a\, bi, a 2 , and b 2 , we 
have determined an estimate of u* and uj 2 from the equations 

E(pv/u) p (ati, 6i;u*(q;i, 61), 0) = 0, (3.5) 
w 2 (a 1 ,6 1 ,a2,6 2 ;p) = E(dl3y/dv) p (a 2 ,b 2 ;u* p (a 1 ,b 1 ),0). (3.6) 

Note that u*(ai, bi) is determined implicitly by Eq. (|3.5| ), which has been solved numerically 
for each value of a\ and b\. 

Then, we have considered sets of approximants such that a^ 2 G [a — Aa,a + Aa] and 
^1,2 € [6 — Ab,b + Ab]. The final estimate was obtained averaging over all integer values 
of «i, a 2 , bi, and b 2 belonging to these intervals. The results for iV = 3 and several 
choices of the parameters are reported in Table [VI]. In order to obtain a final estimate, 
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we should devise a procedure to determine "optimal" values for the parameters a, b, Aa, 
and Ab. Reasonable values of a and b can be obtained by requiring that the estimates are 
approximately independent of the number of terms one is considering. It is less clear how 
to determine the width of the intervals Aa and Ab. Indeed, the results are stable, within 
the quoted errors, for many different choices of these two parameters, while the standard 
deviation of the estimate, which we use as an indication of the error, strongly depends on 
the choice one is making. In order to have reasonable error estimates, we have compared our 
results for u* with the estimates obtained from the analysis of the same series by different 
authors. For N = 3, Guida and Zinn- Justin quote 1.390(4) [pB" |. We have therefore chosen 



our parameters so that we reproduce their errors. More precisely, we choose Aa = 2 and 
Ab = 3, and quote the error in the final results as two standard deviations. With this choice, 
we obtain u* = 1.393(4), which agrees with the previous estimate and has the same error. 
Results for other values of N are reported in Table |V11] . Again, one can verify the good 
agreement of our estimates of u* with the results of Ref. [j6^] . 

In Table |V11] we also report our results for uo 2 . The 0(2)-symmetric point is stable, since 
uo 2 = 0.103(8) > 0. On the other hand, the symmetric point is clearly unstable for N > 4. 
For N = 3, the analysis gives uo 2 = —0.013(6) < 0, so that the fixed point is unstable and 
N c < 3. To better understand the reliability of the results, in Fig. ^|we show the distribution 
of the estimates of uj 2 when we vary a 1; b±, a 2 , and b 2 in the chosen interval. It is evident 
that the quoted error on uj 2 is quite conservative and that the results become increasingly 
stable as the number of loops increases. 

We have then determined N c , defined as the value of N for which uo 2 = at the O(N)- 



symmetric fixed point [73|. The computation was done as before. For each value of the four 



parameters ai, b x , a 2 , and b 2 , we computed an estimate of N c , by requiring 

Q 2 {a 1 ,b 1 , a 2 ,b 2 ;p) = 0. (3.7) 

The distribution of the results is shown in Fig. [3]. It is evident that N c < 3, as of course it 
should be expected from the analysis of u 2 at N = 3. The result is increasingly stable as 
the number of orders that are included increases. We estimate 

N c = 2.91(9) (4 loops), 2.91(3) (5 loops), 2.89(2) (6 loops). (3.8) 

We can thus safely conclude that N c = 2.89(4). 

We have then considered the cubic fixed point and studied its stability. Again we used 
four parameters ai, b\, a 2 , and b 2 . For each choice we first computed an estimate of the 
critical point u*(ai, bi, a 2 , b 2 ), v p (ai, b±, a 2 , b 2 ) solving the equations 

E(f3^/u) p (ai,b 1 ]V^(ax,b 1 ,a 2 ,b 2 ),v*(ai,b 1 ,a 2 ,b 2 )) =0, (3.9) 
E((3 1J /v) p (a 2 ,b 2 ]u* p (a 1 ,b 1 ,a 2 ,b 2 ),v*(a 1 ,b 1 ,a 2 ,b 2 )) =0. (3.10) 

Then we determined the elements of the stability matrix Q from 

Q n = E(d/3u/du)p(a 1 ,bi;u^ ) (a 1 ,b 1 , a 2 , b 2 ), v*(ai,h, a 2 ,b 2 )), (3-H) 

tt 22 = E{d^/dv) p {a 2 ,b 2 \ul{a x ^oi 2 M)^%{^iM^ 2 ,b 2 )), (3.12) 

fii2 = Tt p (ai,b 1 ,a 2 ,b 2 )E(l/u ■ d^/ dv) p ^i(a x , h; w*(ai, b u a 2 , b 2 ),v* p {oti, h, a 2 , b 2 )), (3.13) 

^21 = v*(a x , 6i, a 2 , b 2 ) E(l/v ■ dp u /du) p - 1 (a 2 , h; w*(ai, 6 X , a 2 , b 2 ), v* p {a u b u a 2 , b 2 )). (3.14) 
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'2.70 2.75 2.80 2.85 2.90 2.95 3.00 3.05 3.10 

FIG. 3. Distribution of the estimates of N c . It is determined by requiring uj 2 = at the 
0(./V)-symmetric fixed point. 



We computed the eigenvalues of f2, obtaining estimates of the exponents u\ and oj 2 . As 



before, we determined a and b by requiring the stability of the estimates of u*, v*, u\, 



and 



cu 2 , when varying the order of the series. Errors were computed as before. The results 
are reported in Table [VII]] . They show that for N > 3 the cubic fixed point is stable. 
Note that for iV = 3 we find v* > 0, which agrees with what is expected in the case of 
a stable cubic fixed point, see Fig. |l|. The distribution of the estimates of v*, obtained 
varying the parameters ai,bi,a 2 ,b 2 , and reported in Fig. f|, shows that the result is quite 
stable. In accordance with this scenario, we find ui 2 = 0.010(4) > 0. The quoted error is 
quite conservative, as it can be seen from Fig. ||] in which we report the distribution of the 
estimates of u 2 . Note that the results become increasingly stable as the number of loops 
increases. One may compare our estimate of uo 2 at the cubic fixed point with the four-loop 
uj 2 = 0.0081, which is fully consistent with ours J74] . 



result of Ref. ffi 
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TABLE VIII. Estimates of the cubic fixed point and of the eigenvalues of the stability matrix 



for various values of N and orders p of the perturbative series. Our final results correspond to 
p = f . 



N 




p 


u* 




V* 




L0 2 


3 




4 


1.36(11) 




0.04(13) 


0.780(11) 


0.011(31) 






5 


1.328(26) 




0.089(28) 


0.777(5) 


0.009(6) 






6 


1.321(9) 




0.096(10) 


0.781(2) 


0.010(2) 






f 


1.321(18) 




0.096(20) 


0.781(4) 


0.010(4) 


4 




4 


0.907(72) 




0.606(82) 


0.711(77) 


0.144(78) 






5 


0.883(17) 




0.639(17) 


0.804(48) 


0.049(40) 






6 


0.881(7) 




0.639(7) 


0.781(22) 


0.076(20) 






f 


0.881(14) 




0.639(14) 


0.781(44) 


0.076(40) 


8 




4 


0.448(49) 




1.138(65) 


0.695(85) 


0.211(93) 






5 


0.440(14) 




1.140(16) 


0.831(77) 


0.098(47) 






6 


0.440(6) 




1.136(5) 


0.775(44) 


0.149(33) 






f 


0.440(12) 




1.136(10) 


0.775(88) 


0.149(66) 


oc 




4 


0.182(15) 




1.422(24) 


0.744(95) 


0.185(15) 






5 


0.173(6) 




1.424(7) 


0.783(31) 


0.177(6) 






6 


0.174(3) 




1.417(3) 


0.790(9) 


0.178(3) 






f 


0.174(6) 




1.417(6) 


0.790(18) 


0.178(6) 




TABLE IX. 


Estimates of the critical exponents at the cubic fixed point for various values of TV" 


and order p of the perturbative series. The column "final" reports our final results. 


The results of 


Ref. 46], obtained by a Pade-Borel analysis of the four-loop series, 


are reported in the last column. 


N 






p = 4 


p = 5 


p = 6 


final 


Ref. |46] 


3 


7 




1.399(46) 


1.390(11) 


1.390(6) 


1.390(12) 


1.3775 




V 




0.710(26) 


0.706(7) 


0.706(3) 


0.706(6) 


0.6996 




V 




0.0338(90) 


0.0325(20) 


0.0333(13) 


0.0333(26) 


0.0332 


4 


7 




1.415(41) 


1.403(12) 


1.405(5) 


1.405(10) 


1.4028 




V 




0.718(23) 


0.714(7) 


0.714(4) 


0.714(8) 


0.7131 




V 




0.0323(73) 


0.0305(18) 


0.0316(11) 


0.0316(22) 


0.0327 


8 


7 




1.401(38) 


1.402(12) 


1.404(5) 


1.404(10) 


1.4074 




V 




0.709(21) 


0.711(7) 


0.712(3) 


0.712(6) 


0.7153 




V 




0.0292(61) 


0.0296(18) 


0.0306(10) 


0.0306(20) 


0.0324 


oc 


7 




1.400(21) 


1.395(9) 


1.396(7) 


1.396(14) 






V 




0.709(12) 


0.707(5) 


0.708(4) 


0.708(8) 






V 




0.0287(35) 


0.0294(11) 


0.0304(10) 


0.0304(20) 





17 




FIG. 4. Distribution of v* of the cubic fixed point for N = 3. 



We computed the exponents at the cubic fixed point, using Eqs. ( |2.9| ), ( |2.10| ), and 
( 2. The results are reported in Table fX[ Note that, for N = 3, the exponents at 
the cubic fixed point do not differ appreciably from those of the isotropic model. A recent 
reanalysis of the fixed- dimension expansion of the three-dimensional 0(iV)-symmetric 
models obtained the following estimates for N = 3: 7 = 1.3895(50), v = 0.7073(35), and 
i] = 0.0355(25). These results should be compared with the critical exponents at the cubic 
fixed-point: 7 = 1.390(12), v = 0.706(6), and 77 = 0.0333(26). The difference is smaller than 
the precision of our results, which is about 1% in the case of 7 and v. 

We also checked the exact predictions for the exponents in the large- N limit. For iV — > 00 
with v and u fixed, we have j3y(u,v) = /3(v)\ lsing , so that v* = v*j. Indeed, our estimate of v* 
is in agreement with the estimate obtained with the same method of analysis in Ref. ||50|| : 
v*j = 1.416(5). We have then compared our estimates rj = 0.0304(20), v = 0.708(8), and 
7 = 1.396(14) with the predictions 
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0.3 



0.2 



0.1 



4 loops 

5 loops 

6 loops 



0.0 




-0.04 



-0.02 



0.00 



0.02 



0.04 



0.06 



CO, 



FIG. 5. Distribution of the estimates of 0J2 at the cubic fixed point for N = 3. 



r) = rn = 0.0364(4) 



(3.15) 



7 



1-aj) 

7/ 
[1-aA 



0.7078(3), 
1.3899(7). 



The Ising model results have been taken from Ref. ||. There is a substantial agreement, 
although we note a small discrepancy for rj. This small difference is not unexpected. Indeed, 
the estimate of rj for the Ising model obtained from the fixed-dimension expansion shows a 
systematic discrepancy with respect to high-temperature and Monte Carlo results. 

Finally, we checked the prediction LO2 = —ot-i/vi at the Ising fixed point. It is immediate 
to verify that 



cu 2 



(3.16) 
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is independent of N. Numerically we find 

uj 2 = -0.174(31) (4 loops), co>2 = -0.178(7) (5 loops), uj 2 = -0.177(3) (6 loops). 

(3.17) 

Our final estimate is therefore uo 2 = —0.177(6) which should be compared with the prediction 
uj 2 = -ai/ui = -0.1745(12). 

IV. ANALYSIS OF THE e-EXPANSION FIVE-LOOP SERIES 

In this Section we consider the alternative field-theoretic approach based on the e- 
expansion. The /3-function and the exponents in the cubic model are known to 0(e 5 ) [flUf . 
These series have already been the object of several analyses using different resummation 
methods [^0H43| , [i6l . Ref. |42| studies the stability of the 0(./V)-symmetric and of the cubic 



fixed points for iV = 3. They rewrite the double expansion in u and v in terms of g = u + v 
and of the parameter 5 = v/(u + v). Then, each coefficient in the expansion in terms of 
5, which is a series in g, is resummed using the known large-order behavior. Since v*, and 
hence 5*, is small, one expects the resulting series to be rapidly convergent near the fixed 
point, and therefore no additional resummation is applied ||41]|. For the cubic fixed point, 



Ref. |42| quotes u 2 ~ 0.0049, 0.0085, 0.0021, obtained respectively from the analysis of 
the three-, four-, and five-loop series. Similar results (with the opposite sign) are found for 
the 0(iV)-symmetric fixed point. These estimates are compatible with a stable cubic fixed 
point, but the observed trend towards smaller values of \ui 2 \ leaves open the possibility that 
the estimate of uo 2 will eventually change sign, modifying the conclusions about the stability 



of the fixed points. The same expansion has been analyzed in Ref. [43[ using a Pade-Borel 
resummation technique. The authors report the estimate N c pa 2.86, but again the uncer- 
tainty on this result is not clear. We will present here a new analysis of the five-loop series 
using the method presented in the previous Section and the value of the singularity of the 
Borel transform given below. 

As before, we will consider expansions in u with z = v/u fixed. The singularity closest 
to the origin of the Borel transform is given by 

' - (1 + z) for z > 0, (4.1) 



1 

u b 



1 + ±z) for z < 0. 



The fixed-point values of u and z associated to the 0(iV)-symmetric and to the cubic fixed 
points are respectively 

«* = ^e + 0(e 2 ), z* = 0, (4.2) 

and 

u*=±e + 0(e 2 ), z* = ^ + 0(e). (4.3) 
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TABLE X. e-expansion estimates of the critical exponents u\ and 0J2 at the cubic and at the 
0(A r )-symmetric fixed points for N = 3. The number p indicates the number of loops included 
in the analysis, the columns labelled by a and b indicate the intervals of a and b used, "final" 
indicates our final estimate from the given analysis. 

a b p = 3 p = 4 p = 5 final 



0(iV)-symmetric fixed point 





[-0.5,1.5] 
[-0.5,1.5] 
[-1.5,2.5] 


[11,17] 
[8,20] 
[11,17] 


0.799(20) 
0.804(36) 
0.830(98) 


0.790(8) 

0.784(20) 

0.784(19) 


0.795(4) 

0.801(14) 

0.795(10) 


0.795(8) 

0.801(28) 

0.795(20) 




[0.0,2.0] 
[0.0,2.0] 
[-1.0,3.0] 


[12,18] 
[9,21] 
[12,18] 


0.004(6) 
0.005(9) 
0.022(32) 


0.000(3) 

-0.002(5) 

0.000(7) 


-0.003(2) 
-0.003(2) 
-0.003(4) 


-0.003(4) 
-0.003(4) 
-0.003(8) 


cubic fixed point 




[-0.5,1.5] 
[-0.5,1.5] 
[-1.5,2.5] 


[12,18] 
[9,21] 
[12,18] 


0.789(16) 
0.793(31) 
0.82(10) 


0.797(6) 

0.793(13) 

0.794(20) 


0.796(2) 
0.799(7) 
0.793(6) 


0.796(4) 

0.799(14) 

0.793(12) 


uJ 2 


[-0.5,1.5] 
[-0.5,1.5] 
[-1.5,2.5] 


[11,17] 
[8,20] 
[11,17] 


0.005(12) 
0.004(13) 
-0.004(17) 


0.006(4) 
0.006(4) 
0.002(7) 


0.007(2) 
0.006(2) 
0.006(3) 


0.007(4) 
0.006(4) 
0.006(6) 



TABLE XL Estimates of the critical exponent L02 at the 0(A r )-symmetric fixed point for several 
values of N. The number p indicates the number of loops included in the analysis, the columns 
labelled by a and b indicate the interval of a and b used, "final" reports our final estimate. 



N 


a 


b 


p = 3 


p = 4 


p = 5 


final 


2 


[0.0,2.0] 


[9,21] 


0.116(12) 


0.113(5) 


0.114(2) 


0.114(4) 


3 


[0.0,2.0] 


[9,21] 


0.005(9) 


-0.002(5) 


-0.003(2) 


-0.003(4) 


4 


[0.5,2.5] 


[8,20] 


-0.088(11) 


-0.103(7) 


-0.105(3) 


-0.105(6) 


8 


[1.5,3.5] 


[7,19] 


-0.374(24) 


-0.395(12) 


-0.395(4) 


-0.395(8) 



Since only the leading term in the e-expansion is relevant for the determination of the 
singularity e b p9|,|50|, using Eqs. (|4^) , ( fl~3l) , and ([□]), we find 



N + 8 



(4.4) 



for the 0(AQ-invariant fixed point p9| . ^0| , and 

3JV2 for 
for 



efe = - 



4(JV-1) 
3N 



N < 4, 
> 4, 



(4.5) 



for the cubic fixed point. 

The method of analysis is the one explained in the previous Section with two simplifi- 
cations: first, we consider expansions in one variable only; moreover, the value at which we 



21 



TABLE XII. Estimates of the critical exponents at the cubic fixed point for several values of 
N. The number p indicates the number of loops included in the analysis, the columns labelled by 



a and b indicate the intervals of a 


and b used, "final" 


reports our final estimate. 




N 


a 


b 


p = 3 


p = 4 


p = 5 


final 


UJi 


3 


[-0.5,1.5] 


[9,21] 


0.793(31) 


0.793(13) 


0.799(7) 


0.799(14) 


4 


[-1.0,1.0] 


[10,22] 


0.790(20) 


0.787(8) 


0.790(4) 


0.790(8) 


8 


[-1.0,1.0] 


[8,20] 


0.779(16) 


0.782(5) 


0.786(3) 


0.786(6) 


oo 


[-1.0,1.0] 


[10,22] 


0.772(22) 


0.789(11) 


0.802(9) 


0.802(18) 




3 


[-0.5,1.5] 


[8,20] 


0.004(13) 


0.006(4) 


0.006(2) 


0.006(4) 


4 


[-1-5,0.5] 


[11,23] 


0.073(15) 


0.078(5) 


0.078(2) 


0.078(4) 


8 


[-1.5,0.5] 


[11,23] 


0.154(9) 


0.155(4) 


0.155(1) 


0.155(2) 


oo 


[-2.5,-0.5] 


[11,23] 


0.210(17) 


0.208(6) 


0.202(4) 


0.202(8) 


7 


3 


[1.5,3.5] 


[3,15] 


1.370(29) 


1.375(8) 


1.377(3) 


1.377(6) 


4 


[1.5,3.5] 


[3,15] 


1.414(20) 


1.421(7) 


1.419(3) 


1.419(6) 


8 


[0.5,2.5] 


[3,15] 


1.420(22) 


1.424(8) 


1.422(3) 


1.422(6) 


oo 


[0.0,2.0] 


[3,15] 


1.394(27) 


1.399(11) 


1.399(4) 


1.399(8) 


V 


3 


[1.5,3.5] 


[3,15] 


0.695(16) 


0.699(5) 


0.701(2) 


0.701(4) 


4 


[1.5,3.5] 


[3,15] 


0.717(11) 


0.723(5) 


0.723(2) 


0.723(4) 


8 


[1.0,3.0] 


[3,15] 


0.723(7) 


0.724(3) 


0.723(1) 


0.723(2) 


oo 


[1.0,3.0] 


[3,15] 


0.713(4) 


0.711(2) 


0.711(1) 


0.711(2) 


V 


3 


[1.0,3.0] 


[3,15] 


0.0319(61) 


0.0359(16) 


0.0374(11) 


0.0374(22) 


4 


[0.5,2.5] 


[4,16] 


0.0319(23) 


0.0339(10) 


0.0357(9) 


0.0357(18) 


8 


[0.5,2.5] 


[3,15] 


0.0301(31) 


0.0336(11) 


0.0349(8) 


0.0349(16) 


oo 


[0.5,2.5] 


[3,15] 


0.0296(36) 


0.0332(14) 


0.0349(11) 


0.0349(22) 



should compute the expansion is known (e = 1). Each series is resummed using several differ- 
ent values of a and b belonging to the intervals a G [a — Act, a + Act] and b G \b — Ab, b+Ab]. 
The parameters ct and b were chosen as before, while, after several trials, we fixed Aa = 1 
and Ab = 6. The error in the final results is always two standard deviations. The depen- 
dence of the results on the different choices of ct and b is shown in Table [X[ The final results 
are reported in Tables [XI| and [XII . 

The final results are in reasonable agreement with the estimates of the previous Section. 
However, the instability of the 0(iV)-symmetric fixed point is less clear in the e-expansion: 
indeed we find uj 2 = —0.003(4) at the symmetric fixed point and lj 2 = 0.006(4) at the cubic 
fixed point. This is not surprising: for the 0(iV)-symmetric model the five-loop e-expansion 
gives results that are less precise than the six-loop expansion in fixed dimension. We also 
mention that for iV = oo we obtain an estimate of T], r\ — 0.0349(22), that is in perfect 
agreement with the theoretical prediction 77 = 0.0364(4). 
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Finally we compute N c . We will determine it following two different strategies. A first 
possibility is to consider u 2 and determine the value of N for which the three-dimensional 
estimate of uo 2 vanishes, which is exactly what we did in the fixed-dimension expansion. At 
the symmetric fixed point we find 

N c = 3.07(9) (3 loops), 2.99(6) (4 loops), 2.98(2) (5 loops), (4.6) 

while at the cubic fixed point 

iV c = 3.04(10) (3 loops), 2.95(6) (4 loops), 2.96(3) (5 loops). (4.7) 

Averaging the two results and reporting the error as two standard deviations, we obtain 
N c = 2.97(6), which is in agreement with the analysis of the same series of Ref. [ID 



N c ~ 2.958. It is also consistent with the results of the analysis of the eigenvalues u 2 
for N = 3: indeed we found u 2 > and u 2 < at the cubic and at the symmetric point 
respectively, implying iV c < 3. However, the error bars on uj 2 are such that the opposite 
inequalities are not excluded. Analogously here, although iV c < 3, the error is such that it 
does not exclude N = 3. In practice, from this analysis alone, it is impossible to conclude 
safely that N c < 3. 

A second possibility consists in solving the equation uo 2 = perturbatively, obtaining for 
N c an expansion in powers of e. The result is independent of which fixed point is chosen. 
Unfortunately, the singularity of the Borel transform of N c (e) is not known and therefore, 



we used the Pade-Borel method. The analysis of the series already appears in Ref. [46 
and will not be repeated here. Instead, we will try to make use of the fact that N c = 2 for 
d = 2, performing a constrained analysis. The method has already been applied in many 



instances |]75| , |7^ , |7^ , |77l , |7T| , providing more precise estimates of the critical quantities. The 



method consists in rewriting 

jV c (e) = 2 + (2-e)A(e), (4.8) 

where 

w-W- (4 - 9) 

The quantity A(e) is expanded in powers of e and then resummed. One can verify that 
the coefficients of the expansion of A(e) are uniformly smaller than the coefficients of the 
original series, by a factor of two approximately. Therefore, one expects a corresponding 
improvement in the error estimates. We considered Pade's [2/1], [1/2], [2/2], and [3/1] 
and several different values of the parameter b between the "reasonable" values and 20. 
Pade's [1/2] and [2/2] have a singularity on the real positive axis for 6^7: these cases are 
of course excluded from consideration. At four loops, we find that the approximant [2/1] 
gives estimates varying between 2.82 and 2.87, while the approximant [1/2] is stable giving 
iV c pa 2.82. At five loops, the approximant [3/1] varies between 2.85 and 2.92, while the 
approximant [2/2] is more stable and gives N c pa 2.83. To appreciate the improvement of the 
results due to the constrained analysis, we report the corresponding variation of the estimates 
for the original series: for < b < 20 the approximants [2/1], [1/2], [2/2], and [3/1] give 
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estimates varying in the intervals 2.80 < A c < 2.92, 2.90 < N c < 2.96, 2.90 < N c < 2.95, 
2.86 < N c < 3.03. A conservative final estimate is N c = 2.87(5). This result is lower 
than that of the previous analysis, but still compatible with it. The constrained analysis is 
in much better agreement with the results obtained in the fixed-dimension expansion, and 
clearly supports the claim that iV c < 3. 
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APPENDIX A: CUBIC FIXED POINT IN THE e-EXPANSION 

We report here the e-expansion computation of the four-point renormalized couplings u* 
and v* defined in Sec. [II|at the cubic fixed point, following Refs. |70|,|7TJ. We have calculated 



the three-loop expansion of these two quantities. The final result can be written as 



, v v*(e)=J2v n e\ (Al) 

diV n=0 n=0 



where 



u = 1, (A2) 
_ (A-1)(19A-106) 

27N 2 > ^ 

2 (-11236 + 22540 A- 14181 A 2 + 2446 A 3 + 107 A 4 ) 
U2 ~ 729 A 4 

A (A- 1) (86 + 7 A) 4 (14-7A-6A 2 + 2A 3 ) C(3) 



U 3 



81 A 2 9 A 3 

8 H (10 - 5 A - 2 A 2 ) A (-9116 + 16328 A - 7863 A 2 + 209 A 3 + 37 A 4 ) 



27 A 3 729 A 4 

1 —— f47640640 - 134959200 A + 141439956 A 2 - 64950380 A 3 

157464 A 6 V 

+ 11140557 A 4 - 136770 A 5 + 11821 A 6 ) 
(16-8A-8A 2 + 3A 3 ) vr 4 (A — 1) (86 + 7 A) ( 7b A + Qi) 



405 A 3 162 A 2 

4 (-68 + 54 A - 4 A 2 + 3 A 3 ) Q 2 40 (18 - 5 A - 2 A 2 - 6 A 3 + 2 A 4 ) C(5) 



27 A 3 27 A 4 

-59360 + 93392 A - 27220 A 2 - 14140 A 3 + 5083 A 4 + 346 A 5 ) ((3) 

486 A 5 ' ' ' 



and 
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(N - l)(-424 + 110JV + 17iV 2 ) 
Vl ~ 2W ' (A7) 

25 A (N - 4) (JV - 1) (2 + N) 4 (2 + N) (28 - 32 JV + 6 TV 2 + iV 3 ) C(3) 

V2 ~ 81 iV 3 9N 1 

2 (—44944 + 93764 iV — 61408 N 2 + 10951 iV 3 + 565 iV 4 + 100 A^ 5 ) 
+ 7^ > ( A8 ) 

8H (N-A) (8-AN-N 3 ) 

V3 ~ 27 

A (-21200 + 31540 N - 6536 N 2 - 5953 N 3 + 581 N 4 + 353 N 5 ) 

729 N 5 

+ 1 —— (190562560 - 555117760 N + 598171440 N 2 - 281214884 N 3 

157464 N 7 V 

+47944832 N 4 + 267603 N 5 - 77234 N 6 + 23315 N 7 ) 

(68 -42A^-32A^ 2 + 14A^ 3 + A^ 4 ) vr 4 _ 25 (N-A) (N — 1) (2 + N) A + Q%) 

405 N 4 162 A^ 3 

4 (N - 4) (-24 + 22N -5N 2 + 2N 3 ) Q 2 



9N 4 

40 (72 -30A^-llA^ 2 -18A^ 3 + 7A^ 4 + A^ 5 ) C(5) 



27 

(-237440 + 401536 N - 151152 N 2 - 36648 N 3 + 15662 N 4 + 2277 N 5 + 68 N 6 ) C(3) 

486 A^ 6 ' 

(A9) 



Here "fE is Euler's constant. We have also introduced the following numerical constants [71 



A = 1.171953619344729445... 
Qi = -2.695258053506736953... 
Q 2 = 0.400685634386531428... 

H = -2.155952487340794361... (A10) 



A standard analysis J70| gives: 



N = 3 
N = A 
N = 8 



u* = 1.416(10), v* = -0.03(14) ; (All) 

u* = 0.971(19), v* = 0.58(9) ; (A12) 

u* = 0.455(73), v* = 1.14(9). (A13) 



The estimates of v* are in reasonable agreement with the results of Section |T|, although 
they are much less precise, as it should be expected since here we are analyzing a shorter 
series. On the other hand, the estimates of u* strongly disagree with the quoted error bars. 
The estimate for A^ = 3 is the worst one: indeed, the six-loop fixed- dimension expansion 
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predicts u* = 1.321(18). However, the quoted errors (based, as usual, on the stability of the 
estimates when changing the parameters b and a) seem to be largely underestimated. For 
instance, the three-loop series for u* at the 0(A r )-symmetric fixed point gives u* = 1.39(7) 



71fl : the error is in this case seven times larger! There is no reason to believe the error on 
the estimate of u* at the cubic point to be much smaller than that at the isotropic one. 
Moreover, it is difficult to accept that an expansion truncated at three loops might give a 
more precise result than the six-loop fixed- dimension expansion. Thus, the previous error 
estimates should not be trusted and the observed stability is probably accidental. If we 
assume that the correct error is of order « 0.07 as in the 0(3) case, then all results are in 
agreement. Note that the errors on v* are instead of the expected order of magnitude. 
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